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OSCILLATIONS  OF  A  RESERVOIR  WITH  A  FLUID  ON  THE  FREE  SURFACE  OF  WHICH 
A  MEMBRANE  IS  LOCATED 


L.  V.  Dokuchayev 
(Kalingrad,  Moscow.  Oblast ' ) 


The  creation  of  large-capacity  reservoirs  requires  the  careful 
analysis  of  the  possible  resonance  excitation  of  wave  motions  of  the 
liquid  by  wind  or  seismic  pulsed  loads.  A  number  of  reports  are  con¬ 
cerned  with  this  problem,  e.g.,  [1,  2],  A  rib,  a  special  membrane, 
or  a  thick  plate  that  covers  the  free  surface  of  the  liquid  can  be 
used  to  limit  the  liquid's  mobility.  The  dynamics  of  this  type  of 
system  are  considered  below. 

So  far,  a  specific  shape  of  reservoir  (sphere,  cylinder)  has  been 
studied  in  all  of  the  works  on  the  dynamics  of  shells  with  a  liquid. 
This  report  considers  oscillations  of  a  rigid  reservoir  on  an  elastic 
base.  The  constant  structure  of  the  differential  equations  of  the 
membrane  makes  it  possible  to  obtain  the  dynamic  characteristics  of 
the  "housing  -  liquid  -  membrane”  system  in  analytical  form  for  the 
general  case  of  the  configuration  of  a  reservoir. 

1.  We  will  consider  an  arbitrary  configuration  of  a  reservoir 
with  a  liquid,  on  the  free  surface  Sq  of  which  a  smooth  plate  that 
is  fastened  on  its  periphery  C  is  placed.  We  will  designate  the  re¬ 
mainder  of  the  wet  surface  of  the  reservoir  by  S  (Fig.  1). 


We  will  assume  that  the  liquid  is  perfect  and  incompressible, 
and  we  will  consider  the  case  of  an  isotropic  plate. 

Suppose  that  a  reservoir  filled  with  a  liquid 
rotates  at  an  angular  velocity  fl»  relative  to  the 
elastic  base  and  moves  at  speed  ^  under  the 
action  of  seismic  or  wind  loads.  We  will  desig¬ 
nate  the  displacement  of  the  plate  particles  on 
the  normal  by  w,  and  the  displacement  potential 
of  the  liquid  -  by  i 

The  plate  has  a  constant  flexural  rigidity 

Fig.  1.  of  D  and  is  subjected  to  tensile  stresses  T  in 

its  middle  surface. 

Then,  in  the  case  of  small  oscillations  of  the  "housing  -  plate 
-  liquid"  system,  we  obtain  the  following  linear  system  of  equations 
in  partial  derivatives: 

Ax-0  (1.1);  a*/*'L=(i,  +  dxld  .  I  (1.2; 

**  5 

?m&m&9/dP  +  Dv*v,m  —  T&t»  =  p;  (1.3) 

v*i5t<oo;  •jc«a./*.|c  =  o.  (1.H) 

Equation  (1.1)  follows  from  the  condition  of  the  incompressibility 
of  the  liquid.  Equation  (1.2)  consists  of  the  conditions  under  which 
the  liquid  does  not  flow  through  the  surface.  Equation  (1.3)  is  the 
equation  of  the  oscillations  of  the  plate;  (1.4)  -  the  conditions  of 
the  boundedness  of  the  oscillations  and  the  boundary  conditions  of 
the  rigid  attachment  of  the  plate;  v*— opr  the  external  normal  to  the 
surface  of  a  liquid  occupying  volurtie  Q;  n  —  opt  the  external  normal 
to  the  contour  C;  the  values  and  designate  the  density  and 

thickness  of  the  plate,  respectively.  The  transverse  load  p  to  which 
the  plate  is  subjected  from  the  side  of  the  liquid  can  be  determined 
with  precision  down  to  the  constant  by  using  the  linearized  Lagrange- 
Cauchy  integral  according  to  the  formula 

#  — — f  — g  (*i  +  +  *••*)!.  (1.5) 
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J,1  JH  .'i.i.'jr  rj  i.t  r  «■.  jvpj  a  rj .-.  if.  *.  -’.  -r  r.  y.  r 


Ktt 


where  t  -  the  density  of  the  liquid;  g  -  the  acceleration  of  gravity. 

We  will  consider  the  homogeneous  problem  of  free  oscillations 
of  the  "flexible  plate  -  liquid"  system  when  the  reservoir  housing 
is  immobile,  and  we  will  use  the  solution  of  the  following  boundary 
problem  on  free  oscillations  of  a  liquid  in  a  mass  force  field: 


Af— 0;  i  f /d  V  |3  =  0;  df/d  v  |s#  =  xf. 


(1.6) 


This  problem  has  an  infinite  discrete  range  of  eigenvalues 
*«  and  the  eigenfunctions  ?«  corresponding  to  them. 


According  to  report  [3].  we  will  isolate  the  time  coordinate, 
and  we  will  expand  the  unknown  function  of  into  a  series  according 
to  this  system  of  functions: 


Ni 


3. 


5# 


(1.7) 


We  obtain  the  function  v  from  the  solution  of  the  homogeneous 
boundary  problem  of  elastic  oscillations  of  a  plate,  which,  according 
to  ( 1 . 3)-( 1 • 4) ,  has  the  form: 


L  (o)  =  D  v*  V*  *  —  T  V*  v  +  (p  g  —  p„  dB  «•*)  o  —  ?  •*  2  f  »S 


v®ls, <»;  =««’/*'» Ic-®* 


( i  .a) 


The  general  solution  of  boundary  problem  (1.8)  is  represented 
in  the  form  of  the  linear  combination  of  the  unknown  solutions  p®  of 
the  homogeneous  equation  L(v)  =  0  and  the  partial  solution  *  of 
heterogeneous  equation  (1.8)  in  the  form  of  the  expansion  of  problem 
(1.6)  into  eigenfunctions. 

Using  (1.4),  and  also  the  second  formula  in  (1.7),  we  obtain  the 
solution  to  problem  (1.8)  in  the  form  of  functions  that  form  a  complete 
orthogonal  system: 
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*  ■  m 

+  2  B*k*Vnl9*> 


(1.9) 


where  the  eigenvalues  are  the  roots  of  the  characteristic  equation 


&  *■  (  j£i  *l~a* 


The  following  designations  are  used  here: 


(1.10) 


~ Ic :  Lr=dt>V*n  !c;  =  *•!  “*«**>,, /a  «Ic: 

-iiilV  .  go  ^  J_  [J>**3dc.  q  m  K-B?« 

‘  ff-  *  “  Nl  { s=r-: 


(1.11) 


Using  the  eigenfunctions  v^,  we  can  construct  the  solution  to 
the  system  of  equations  in  partial  derivatives  (1.1),  (1.3)  and  for 
the  case  of  an  oscillating  reservoir.  We  will  represent  the  displace¬ 
ment  potential  of  the  liquid  in  the  form  of  the  sum  of  terms,  the  first 
of  which  corresponds  to  the  displacements  of  the  liquid  in  an  absolute¬ 
ly  rigid  shell,  while  the  last  corresponds  to  the  elastic  deformations 
of  the  plate 


•  •  • 

xt  *t  w + 2  v*  ***’  **•  *• )•/  (o  +  T  f«  (*i.  X),  *,)  n„ 


(1.12) 


Here  the  functions  are  the  Zhukovskiy  potentials  that  satisfy 
the  following  boundary  problem: 

functions  f«  -  the  eigenfunctions  of  boundary  problem  (1.6).  Expansion 
coefficients  (1.12)  are  determined  by  the  following  formula: 


H  _JL  f  wtSs 

"  Nl  i  •' 


We  will  represent  the  unknown  function  w  in  the  form  of  an  in¬ 
finite  expansion  by  types  of  free  oscillations  of  the  plate  (1.9): 


Substuting  expansion  (1.13)  in  equation  (1.3),  which  makes  it 
possible  to  perform  differentiation  by  terms  and  the  integration  of 
the  series,  and  also  considering  the  properties  of  orthogonality,  we 
multiply  (1.3)  by  v^dS  and  integrate  with  respect  to  So-  As  a  result, 

we  obtain  the  infinite  system  of  equations 

s  a 

(**  (?*+  •*  Pt)  +  2  +  +  2  [  (  —  1/  a  ®/]  =0. 

/—i  /-? 

k  =  1,  2,  ...;  the  dots  designate  differentiation  through  time.  The 
square  of  the  plate  oscillation  frequency  wh  is  the  k-th  eigen  number 
of  the  characteristic  equation  (1.10),  and  the  remaining  coefficients 
are  determined  by  the  following  formulae: 


(**  ■=?  2  a*/  *»  =  2 

M  m 

1  *«•  **<*$= 2  ■  5  2  Bn* 

*.a-p  f  «■/  y  Bmtf  f  * 

5.  a-1  «+**  mmi 


(1.14) 


where  the  coefficients 


J  *-17  iS-  *;•-*£*/  -57 

J  «•  (* xT)/ 

■»  S+5# 


tfS 


(1.15) 


are  determined  by  the  methods  described  in  report  [4],  for  example. 
The  coefficients  Bnk  are  calculated  according  to  formula  (1.11). 


The  forces  and  moments  acting  on  the  body  from  the  side  of  the 
liquid  and  the  elastic  plate  are  determined  by  the  pressure  p,  which 
calculated  using  formula  (1.5)  with  consideration  of  (1.12).  Then 
the  equation  of  motion  of  the  "reservoir  -  liquid  -  plate"  system 
can  be  written  as  follows,  with  the  assumption  that  the  axes  of  the 
introduced  coordinate  system  are  the  primary  ones, 


(*•+*)  »t  +  (  -  1)'  (L»+  L)  +  2  */*>*  =  /»/.*  =  2.3; 

A-l 

m 

(/?  +  /,)  '§,  -  (- 1  )‘  (L*  +  L  )  V<  -  B  (!•  +  L)  0,  +  V  [  *o  «  P*  + 

«• 

+  l(  —  *P*  ]  (  ^?  +  /i)®i  +  ^  ^H*P*=,^**» 

*=i 

where  x4*  and  x.u  coincide  with  coefficients  (1.14);  m0  +  ra,  L.0  +  L, 

I°i  +  Ii  -  the  mass,  static  moment  and  inertial  moment  of  the  system 
with  an  absolutely  rigid  membrane;  and  -  the  projections  of  the 
actions  of  the  elastic  base  of  the  reservoir  and  the  external  loads. 

2.  We  can  construct  an  approximate  solution  for  a  reservoir 
with  an  arbitrary  rotational  contour  for  the  boundary  problem  (1.1)— 
(1.4)  from  the  known  range  of  eigenfrequencies  of  problem  (1.6). 

This  solution  precisely  satisfies  the  Laplace  equation  and  the  boundary 
condition  on  the  free  surface  Sq •  It  satisfies  the  boundary  condition 
on  the  wet  surface  S  on  the  average,  which  is  described  in  rather 
great  detail  in  report  [4], 


In  a  cylindrical  coordinate  system,  it  has  the  following  form; 

»/>*  =  ipm  tpm  (x.  r)  =  Xpm  (x)  Ypm  (r)  (p,  m  =  1 .2  ...)  , 


where 


t'pm  <*— *•)/'•  •  Y  =*  ^ m  r L r . 

Vpiid 


Cpm  -  V"»  r»- 


Here  xq  -  the  coordinate  of  the  free  surface  of  the  liquid;  x,. 

-  the  predetermined  range  of  eigenvalues  of  the  boundary  problem  (1.6); 
ro  -  the  radius  of  the  contour  of  the  free  surface  C. 

Instead  of  the  single-parameter  family  of  subscripts  n,  here  we 
will  introduce  the  two-parameter  family  of  subscripts  p  and  m  corres¬ 
ponding  to  the  radial  and  circular  harmonics. 

We  will  consider  an  inertialess  plate  in  the  absence  of  a  mass 
force  field.  In  this  case,  boundary  problem  (1.8)  for  a  circular 


flexible  plate  is  written  as  follows,  with  the  isolation  of  the  cir¬ 
cular  coordinate: 


j  I'.  B„m  Y,m  (r)/^]; 

P“*  I 

lr»r,  —  0;  dom/ dr  | r=r%  —  0, 

Where  ^r^/D;  U  =  {d/rdr)  {rd/d^-m'lr*.  (2.2) 

The  general  solution  to  this  equation  is  the  function 

=  C,  (r/rt)n  +  C,  Zu  (r)  +  -J,  ^  lBpm  ^pm  (')/apmJ. 

where 

7*  (')  =  (Y  '/'•»)//*  (Y):  apm=D'lm['lm  +  Y*  )  P^- 

According  to  formulae  (1.11),  the  coefficients  and  Bpmk 
assume  the  form 


(2.1) 


_o  ^  (w  Ypw) 

^  4.-*+C 


:  b5»,= 


2  -p«  (  ^m  ~  ygm) 


Bpmk  =  (  ^/n  B1  pm  ~  m  pm)  0pm/(  ®m*  ~  apm)  (  ~  n)  ’ 


(2.3) 


where 


Ypm  —^pm  (C pm)/Jm  (5pm):  ”=  Y  l/m  (Y)  /  An  (Y)- 


Considering  the  property  of  orthogonality,  we  come  to  the  con¬ 
clusion  that  only  harmonics  with  the  subscript  m  =  1  can  be  excited 
during  any  movement  of  the  reservoir. 

All  of  the  coefficients  (1.15)  for  the  different  shapes  of  re¬ 
servoirs  encountered  in  practice  are  determined  in  [4].  Therefore, 
in  order  to  create  a  dynamic  diagram  of  the  "reservoir  -  liquid  - 
plate"  system,  it  suffices  to  calculate  the  coefficients  (2.3)  using 
the  Bessel  function,  and  to  find  the  root  of  characteristic  equation 
(1.10). 

3*  As  an  example,  we  will  consider  a  circular  cylinder  with 
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*>  v  v  v  /.v.v.  w. 


•  Vi**  v*  « 


V  V  V 


4 


si*  i 


fcW 


a  unit  radius  ro  =  1  •  Then  the  eigen  numbers  and  eigenfunctions  of 
boundary  problem  (1.6)  will  be  written  as: 

Pi  "=  ipm-  *pm  ~  Epm  ^  *)• 

(t pm  r)Vm  dpm).  (3-D 

where  -  the  root  p  of  equation  Jm  (6)  = 0 


b 

< 

T 

7 

_ _ 

□ 

2 

J  *  n 

Fig.  2. 


The  characteristic  equation  assumes  the 


form: 


~  2  [*/•  'i  *»1  -  «,«']  =  o , 


p-i 


(3.2) 


where  the  following  dimensionless  values  arr  n;.ro 
duced,  with  consideration  of  (2.3): 


Vn  (&.+Y1):  dpm=2rnyZml\  E^,4-  y*i(  l%n~  "^)  Zm+i- 


Bpm*  ■=  dpm  apm  /  (  »>Lk  “  “*■  )•  ~  tm+l  (>')  /  (Y)i  *VP  'o  ~  1  • 


(3.3) 


The  dynamic  characteristics  for  the  cases  f=i  assume  the  form 
(1.14),  where 


•  *  ( ^  **)  »• 


(3.4) 


In  the  limiting  case  of  an  absolutely  elastic  plate,  *-o.  i.e. 
when  there  are  no  forces  in  the  middle  surface  (T  =  0),  formulae 
(3-2)— (3*4)  remain  valid  if  we  use  the  limiting  transition  for  the 
Bessel  function  of  an  imaginary  argument  yZ"»/z„+,  -  2 m(m+i)  in  them. 
Formulae  (3-2)— (3-4)  also  remain  valid  for  an  absolutely  flexible 
plate,  if  we  use  the  limiting  transition  as  Y-^= 

•i*  /  Y*  -  oJu  :  *mly*m+ 1"*’1’ 

Now  we  will  consider  the  case  of  an  absolutely  flexible  inertia¬ 
less  plate  (D  =0)  in  a  field  of  mass  forces  with  gradient  g. 


In  this  case,  equation  (1.8),  written  in  dimensionless  form  for 


a  cylindrical  cavity,  analogously  with  (2.1)  and  (3-2),  can  be  repre¬ 
sented  as  follows: 


U  (*«)  -  2  *  =  m  'i/r-  (3.5) 

P=1 


With  precision  down  to  normalization,  the  general  solution  to 
equation  (3*5)  with  the  boundary  condition  vm( 1 )  =  0  is  the  function 


cmt  (r>  =  -«&  (')  +  ^  — 


(0) 


(3-6) 


Here  a,„ =o*.»  is  the  k-th  root  of  characteristic  equation  (3-2), 
where 


dpm  — 


pm* 


pm* 


(3.7) 


and  the  dynamic  characteristics  will  be  determined  using  formulae 
(1.14),  (3.6) . 


When  there  are  no  mass  forces  (6-0)  at  the  limit,  equations 
(3.7)  are  transformed  into  (3-3)  as  v-^- 

On  the  other  hand,  when  the  mass  forces  predominate  over  the 
surface  tension  forces  (*-oc).  the  coefficient  and  the  square 

of  the  frequency  and  coef f icients  (  1 .  1 4 )  approach  the  following  val¬ 
ues  (3-4): 

^*/®****im’  ^*-*^*i:  W  p**'  (3-6) 


Infinite  series  (3-2)  and  (3.5)  rapidly  converge.  Figure  2  shows 
the  values  of  the  first  four  roots  (k  =  1,  2,  3»  4)  of  equation  (3-2) 
at  m  =  1  and  v-o  for  a  different  number  of  terms  of  the  sum  for 
subscript  p.  It  is  evident  that  the  value  of  the  first  root 
already  converged  with  three  terms  (n  =  3);  it  suffices  to  take  the 
first  five  terms  of  the  sum  (n  =  5)  in  order  to  obtain  the  values  of 
the  first  three  roots  . 


■ 

mm 

B 

mum 

■ 

mm 

wim 

■ 

mm 

B 

ufsm 

mmmm 

u 

warn 

2  l9r 

Fig.  3- 


1  7  3  m 

Fig.  M . 


o  i 

Fig.  5. 


Figure  3  shows  the  dependence  of  the  square  of  the  dimensionless 
frequency  on  the  dimensionless  parameter  y.  determined  by  formula 
(2.2)  for  different  k,  using  a  logarithmic  scale.  The  horizontal 
broken  lines  designate  the  asymptotes  which  the  curves  approach  as 
0;  this  is  the  case  of  the  absence  of  tension  in  the  middle  sur¬ 
face  of  an  elastic  membrane  wil.i  flexural  rigidity.  The  sloping 
broken  lines  represent  the  asymptotes  that  the  curves  approach  as 
y ■*°c.  i.e.,  the  flexible  membrane  is  not  resistant  to  bending,  and 

it  is  subjected  to  tensile  forces.  The  dependence  of  the  square  of 
the  frequency  on  the  value  of  the  number  of  waves  in  the  peripheral 
direction  m  for  different  k  and  y  is  shown  in  Fig.  4. 


The  first  three  figures  show  the  case  of  the  absence  of  a  mass 
force  field.  Figure  5  shows  the  effect  of  its  gradient,  determined 
by  the  dimensionless  parameter  6  using  formula  (3-5).  The  broken 
lines  correspond  to  the  asymptotic  behavior  of  the  curves  at  6  ap¬ 
proaching  zero  and  infinity. 
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